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Abstract 


The paper is devoted to the theoretical investigation of the circular dichroism 
(CD) and optical rotatory dispersion (ORD) spectra of turbid layers with 
large particles of arbitrary shapes. The cases of oriented cylinders and 
spheroids are considered in detail. It is assumed that the optical thickness of 
layers is low and the refractive index of scatterers is close to the refractive 
index of a host medium, which is usually the case for bio-liquids. 

It is found that the simple correction function allows to transform CD 
and ORD spectra measured for turbid layers to that of particles. This result 
could be of importance for studies of living cells molecular structure 


changes during their life cycles. 
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1. Introduction 


Circular dichroism (CD) and optical rotation dispersion (ORD) 
spectra are sources of information on asymmetry of molecular 
structures. For instance, this asymmetry is an important feature 
of living cells. The CD and ORD spectra are used in organic 
chemistry as a routine tool in stereochemical analysis. 

In most cases chiroptical spectroscopy (CS) is applied 
to substances in solution. However, biological media are 
usually highly inhomogeneous due to the presence of various 
macroscopic particles (e.g., red blood cells). Thus, there is 
a task of generalization of the CS on the case of particulate 
media. Some results in this direction have already been 
reported [1, 2]. 

The major result is that CD and ORD spectra of 
substances inside small particles may differ substantially from 
spectra measured for dispersed layers. They depend on the 
refractive indices of particles, their size (and, possibly, on 
the concentration of particles, if the optical thickness of a 
sample is high). The case of chiral spheres can be treated in 
the framework of exact electromagnetic scattering theory [1]. 
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However, most biological particles are nonspherical in 
shape. So, the generalization should be done, allowing for 
nonspherical shapes of particles in an optically active turbid 
layer. This can be achieved in two ways. The first is to 
use known solutions of electromagnetic scattering problems 
for chiral cylinders and spheroids [3,4]. The second is to 
apply the approximate scattering theories (e.g., the Rayleigh, 
Rayleigh—Gans, discrete dipoles and anomalous diffraction 
approximations [1,2,5]). Here we will use the second 
approach, namely the anomalous diffraction theory (ADT) of 
van de Hulst [5]. 

This is mostly due to the following reasons. Firstly, the 
ADT is perfectly suitable for the solution of applied problems 
due to its intrinsic mathematical simplicity. Secondly, the 
range of shapes which can be considered is wide. This is 
not the case for the exact electromagnetic scattering solutions. 
In fact, solutions of the wave equation for most particle 
shapes occurring in nature has not been derived due to the 
mathematical complexity involved. The important thing is that 
the accuracy of approximate equations can be far better than 
the accuracy of experiments performed. Clearly, it depends on 
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the correspondence of particle parameters in the experiment to 
the range of applicability of the approximate theory. 

For instance, it was obtained [1] that an additional 
multiplier is needed to account for the light scattering effects in 
the chiroptical spectroscopy of disperse media with particles, 
which are much smaller than the wavelength à of the 
incident light. The case of arbitrarily shaped particles at 
kd|m — 1| «1, |m— 1| « 1 (d is the maximal size of 
particles, k = 27/24, m = (m; +mp)/(2),mz and mp are 
refractive indices for left-handed and right-handed circularly 
polarized waves correspondingly) was studied by Bohren [6]. 
It was found that light scattering effects can be neglected in 
this particular case. Correspondingly, CD and ORD spectra 
of substances inside such particles coincide with spectra 
measured for dispersed layers (in the framework of the single 
scattering approximation, of course). This allows us to use 
the CS, developed for homogeneous media, in the case of 
particulate media. 

Clearly, the CD and ORD spectra of turbid layers are not 
equal to CD and ORD spectra of particles substances if the 
condition kd|m — 1| « 1 is not satisfied (either due to the 
small value of the wavelength A or large values of d). The aim 
of this paper is to consider this discrepancy for the case of large 
(kd > 1) optically soft (|m — 1| < 1) nonspherical particles, 
where the van de Hulst approximation [5] can be used. 


2. CD and ORD spectra of turbid layers with 
nonspherical particles 


Let us assume that the optical thickness of a medium in 
question is low, which allows us to neglect the contribution 
of multiply scattered photons. Then the disperse medium 
can be considered (at least as far as the light transmission is 
concerned) as an effectively homogeneous medium with the 
‘effective’ refractive index [1,5]: 


M = 1 — 2riNk™° S(0) (1) 


where N is a number of particles in a unit volume and S(0) 
is the so-called amplitude scattering function in the forward 
direction. Clearly, this equation can be applied only at small 
values of |M — 1|. Equation (1) is easily generalized for the 
case of disperse media with optically active particles. Namely, 
two refractive indices are introduced in this case: M@,—for left- 
handed circularly polarized waves and M,—for right-handed 
circularly polarized waves. The values of Mı and M3 coincide 
for isotropic nonchiral particles. They generally differ for 
optically active particles. This difference is usually small 
(|M; — M| « 1). Nevertheless, it is detectable and can 
be used to study molecular asymmetries of substances inside 
small light scattering globules. Thus, it follows [1]: 


Mj (a) = 1 — 2riNk ° SiO) (2) 
and by definition [1] 


g(a) = “Re (My (A) — M3(A)) (3) 


9(A) = “Im (My (A) — M3(A)) (4) 


where (à) and 0 (à) are ORD and CD spectra per unit length 
of a sample. 
Let us introduce the complex number: 


WA) = g(a) — i00). (5) 


Then it follows from equations (2)—(5): 


ia N 
WA) = Gr ASO) (6) 
where 
AS(A) = & (å) — Sı (2). (7) 


Thus, the problem of the calculation of CD and ORD spectra 
of a turbid layer is reduced to the calculation of the difference 
AS(à). The value of A S (å) is small due to the small difference 
Am(A) = my (A) — mR(A). It follows due to this fact: 


0s AS(A) 
S'A) = — = : 8 
A) OM | nom  Am(À) (8) 
Substituting equation (8) into (6), one obtains: 
WA) = SS) WOR) (9) 
~ BV ° 
where V is the average volume of particles, 
Wo(A) = pÀ) — ioà) (10) 
UC 
po) = 7 Re (mL) — mr(A)) (11) 
UC 
(A) = FE (mL) — mr(A)) (12) 
and c = NV is the volumetric concentration of particles. 


One can see from equation (9) that initial spectra Yo(A) = 
(0 (A) —i09(A) indeed change and this change is due to spectral 
variation of the function 4? S’ (A). 

Thus, the main task now is to calculate the dependence 
S(A) and then to consider the derivative S’(4). For this, 
we need to make some assumptions. In particular, we will 
assume, as it is usually done in the framework of the anomalous 
diffraction approximation [5], that particles are optically soft 
(jm — 1| « 1) and large (kd > 1). Then the function S (à) can 
be found in the framework of the van de Hulst (or anomalous 
diffraction) approximation [5]: 


2 
50) == f d-aia a3) 
27 yr 

where & is the projection of a particle on the plane 
perpendicular to the incident light and o (p) = kl(m — 1)&(p) 
is the phase shift of a ray on the length /E(p) where the 
parameter l = V/ © is introduced for the further convenience. 
Equation (13) can be applied only to the case of identical 
particles with a fixed orientation. For more complex cases, 
one should use values (S(A)) averaged over the particle’s 
dimensions and orientations. It follows from equation (13) 
for the derivative 


ik3] 
SQ) => / &(p) exp[ —ik(m — lE(py|dp. (14) 
T JE 
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Thus, one can obtain from equations (9), (14): 


WA) = H(A) Wola) (5) 


where 


HO) = f E0 exp ipto) aq (16) 
s 
where q = p/X, p = kl(m — 1), and S is the unit area. The 
shape of this unit area coincides with the shape of the particle 
projection. It follows from equation (15) 
Y) 


Va) = —— 


n (17) 


Thus, the knowledge of the function H (à) allows us to find 
the intrinsic spectrum (A) of particles from the measured 
spectrum W(A). It is known, however, that H(A) — 1 as 
p — 0 [6]. This means that 


t= f e@ag =", (18) 
This equality holds for any particle shape. 

The function H(A) for a given particle size, shape and 
orientation can be calculated numerically or analytically in 
some cases (see below). For numerical integration, one should 
find the length L(q) = /&(q) of a light ray inside the particle. 
After that the integration (16) can be performed for a particle 
of an arbitrary shape. As a matter of fact, a three-dimensional 
particle is replaced by a two-dimensional amplitude-phase 
screen in the framework of the van de Hulst approximation [5]. 
Thus, the main task is to find the phase shift distribution on 
that screen. 


3. Cylinders 


Let us now consider the case of cylindrical particles. We 
will assume that tops of cylinders are oriented normally to the 
incident light beam. The transverse cross section of a cylinder 
can be arbitrary. Clearly, for the illumination of the cylinder 
along the symmetry axis all phase shifts are equal and the 
equivalent amplitude-phase screen has a constant value of the 
phase shift. Then it follows: € = 1 and 


p=kL(m—1) (19) 


where L is the length of the cylinder. We obtain from 
equation (16): 
H(A) = exp(—ip(a)) (20) 


where we accounted for the fact that §(¢) = 1 in this case. It 
follows from equation (15) 

WA) = YA) exp(—ip(A)). (21) 
Then we have from equation (21) (in correspondence with [6]) 
as p(A) > 0: 


WA) = WA). (22) 


However, at larger values of o(A) considerable differences 
between spectra Y (à) and Wo(A) occur. 
In particular, it follows (see equation (21)): 
P(A) = Go(A) cos p(A) — Oo) sin p (à) (23) 
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O(A) = po) sin o (å) + A9(A) cos p(A) (24) 


if particles are not absorbing. It follows at p(A) = 27n 


(n = 0,1, 2,3,...): 


pA) = gol) O(A) = 0). (25) 
However, one obtains at p(A) = a(n + 3) (n = 0, 2, 4, 6, 8): 
g(a) = —A(A) O(A) = Ø). (26) 


Equations (23), (24) allow to consider the distortion of spectra 
for arbitrary phase shifts p (à). 


4. Ellipsoids 


Let us consider now the case of oriented ellipsoidal particles. 
The function €(q) for this case was found by Lopatin and 
Sidko [7]. It has the following simple form [7]: 


TORE NET] 


It follows from equations (16), (27): 


1 
no=;f afl 7exp( ipy1 7) aa (28) 


(27) 


where 


p = kd(m — 1) d= 3. (29) 


For spheres the value of d coincides with the diameter. It 
follows from equation (28) after integration: 


6i(1—e-”?) 6e? 


p? p? 


3ie™’? 
H(A) = . (30) 


We obtain from equation (30) as p > 0: 
H(A) = (1 — +p — 4p’). 


One can see again that H — 1 as p — 0 as it should be [6]. 
The value of H(A) in equation (30) can be presented in 
the following form 


A(A) = u(à) — iv(A) (31) 
where aa č Be 
sin p cos p sin p 
u(a) = + = 2 (32) 
pP p pP 
3cosp O6sinp 6(1—cosp) 
va) == +t — -~ (33) 
p pr pP 


and we assume that particles are nonabsorbing. Then we have: 
PÀ) = u(A)po(A) — và) A) 


0A) = vA) + u(A)O0(A) 


(34) 


(35) 


which is similar to equations (23), (24). Generally speaking, 
linear relationships (34), (35) hold for any particle shape. 
However, the specific dependencies u(A) and v(A) differ 
depending on the specific shape of particles. We present 
functions u(A) and v(A) for cylinders and ellipsoids in figures 1 
and 2. One can see that multipliers H(A) for cylinders and 
ellipsoids differ considerably. In particular, it follows that 
u —> 0, v > Oas p > œ for ellipsoids. This is not the 
case for oriented cylinders. In practice, of course, there is 
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Figure 1. The dependence of the function u on the phase shift p for 
cylinders (solid curve) and ellipsoids (broken curve). 
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Figure 2. The same as in figure 1, but for the function v(p). 


some absorption of light inside particles. This will lead to the 
damping of oscillations both for cylinders and spheroids. 

For illustrative purposes we present spectra g(A) and 6(A), 
obtained with equations (32)—-(35) for different values of the 
parameter b = 22d(m — 1) in figures 3 and 4. Note, that 
we neglected the imaginary part of m in calculations of the 
parameter b = Ap. Also we neglected possible variations 
of b with A. Spectra An(A) = Re(Am(A)) and Ax (à) = 
Re (Am(A)), which close to those of the poly-L-glutamic acid, 
were taken from [8]. Then we used equations (11), (12) to find 
go(A) and 9(A). The concentration of particles c was assumed 
to be equal 0.001. Finally, we expressed go(A) and 6 (A) in 
degrees, multiplying (11) and (12) by 180/z. The thickness of 
a dispersed layer was assumed to be equal 1 mm. We see that 
the parameter b, which is proportional to the phase shift of a 
ray on the maximal dimension of a particle has a paramount 
importance for spectra g(A) and @(A). At small values of b 
spectra of dispersed layers are very close to those of molecular 
solution of poly-L-glutamic acid. This is in a full accordance 
with data presented in [6]. However, they differ considerably 
from initial spectra go(A) and @9(A) for large values of b. 
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Figure 3. Spectra g(A) at b = 0, 20, 100, 500 and 2500 nm. 


0.10 
0.08 
0.06 
ois fi 
0.02 L 

0.00 È 


-0.02 


0, degrees 


-0.04 h 
-0.06 
-0.08 


-0.10 


180 190 200 210 220 230 240 
à, nanometers 


Figure 4. The same as in figure 3, but for spectra 0 (À). 


Generally, we have for large values of b > 00 : g(A) > 0 
and @(A) — 0. This means that the information, contained 
in initial spectra is lost due to the light scattering process. 
However, at small and intermediate values of b the situation 
is not so hopeless (see figures 3 and 4) and information on 
the arrangement of chemical groups for substances condensed 
inside small particles, which is responsible for the particular 
shapes of spectra @o(A) and 6(A), can be retrieved in principle. 
For this, one needs to have an additional information on the 
shape, size and internal structure of particles under study. 


5. Conclusion 


The influence of the shape of particles on the circular dichroism 
and optical rotatory dispersion spectra of turbid layers was 
studied. The investigation was performed in the framework of 
the van de Hulst approximation [5], which is valid for optically 
soft (jm — 1| « 1) large (ka > 1) particles. It was found that 
intrinsic spectra Yo (A) are related to measured spectra Y (A) by 
equation (15), where the function H (A) is given by integral (16) 
for particles of any shapes. The function &(q) in equation (16) 
depends on the shape of particles. It is equal to 1 and 
1.5/1 = q for oriented cylinders and ellipsoids respectively. 
The dependence &(q) for irregularly shaped particles does not 
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[2] 


have such a simple analytical form. However, this function 


can be calculated numerically due to the fact that the product 
1E(q) has a simple geometrical meaning. This is the length of 
a ray inside of a particle. 

A major conclusion is that intrinsic CD 69(A) and ORD 
go(A) spectra of substances bounded inside optically soft 
large particles differ considerably from those measured for 
turbid layers. This difference disappears only for small phase 
shifts [6]. 

The results presented in this paper can be used for 
the development of chiroptical spectroscopy of biological 
media [9-12]. 
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